Moire patterns on STM images of graphite from surface and subsurface rotated layer 
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We have observed with STM moire patterns corresponding to the rotation of one graphene layer on 
HOPG surfaces. The moire pattern domains were characterized by rotation angle and extension in 
the plane. Additionally, by identifying the border domain and defects we can discriminate between 
moire patterns due to rotation on the surface or subsurface layers. For a better understanding of 
moire pattern formation we have studied, from a theoretical point of view, an array of three graphene 
layers where the top or middle layer appears rotated around the stacking axis. By applying first- 
principles total-energy electronic structure calculations to this array we obtain the spatial charge 
density and its corresponding STM image. We compare the experimental and theoretical results 
and show the strong influence of rotations both in surface and subsurface layers for moire pattern 
formations in HOPG STM images. 



PACS numbers: 68.37.Ef, 71. 15. Mb, 73.22.Pr 

The carbon structures have attracted a lot of inter- 
est from several groups. At first, a few decades ago and 
from the fundamental point of view, the STM observa- 
tion of the HOPG surface with atomic resolution was 
reached even at under different experimental conditions 
[IHSj. Notwithstanding, a conclusive explanation for the 
observed triangular and/or honeycomb surface lattice has 
yet to be reached [IHZ]- In the nineties the research on 
carbon structures shifted the attention to the superlat- 
tices [SUTO] , and recently these fields were reactivated by 
the possibility of getting a single and multi graphene layer 
islands onto different substrates such as: Palladium [11 , 
Iridium [12], Ruthenium [13], Copper [I4, or Platinum 

All the hypotheses to explain the superstructures have 
to do with rotational moire patterns resulting from the 
overlap between a missoriented top graphene layer and 
the underlying single graphite crystal. This idea was first 
proposed by Kuwabara et al. [8 , but it became stronger 
when Xhie et al. [10 showed experimental evidence in 
its favour and suggested an explanation based on the 
theory of STM image formation for HOPG surfaces [16]. 
According to this theory only carbon atoms located in /3- 
type sites (these atoms have no neighbors in the adjacent 
layers) are visible for the STM, because they present a 
local density of states (LDOS) larger than those atoms 
located in the a-type sites (which have neighbors in the 
layers above and below). A rotation in the top layer then 
produces regions with different concentrations of atoms 
in the /3-type sites, and those regions where such concen- 
trations are larger will appear brighter at the STM tip. 
However, some points still remain controversial. Because 
it is widely accepted that STM can only render images 
of the surface of a solid, consequently, a relative rotation 
between the top layer and the underlying graphite could 
not be shown directly in the corresponding STM image. 
Then again, there is an alternative explanation for the 



superstructures [9 , and although it is also based on the 
moire patterns, it considers results from first-principles 
calculations [17 to show that current density maxima are 
over regions with a high concentration of atoms in the a- 
type sites. Recently, and in favour of this explanation, it 
was shown from ab initio and tight-binding combined cal- 
culation [18], that for small rotation angle {0 < 3°), the 
electrons with energy close to Dirac point are localized 
in a-type sites. 

In this paper we present experimental evidence of the 
observation of moire patterns induced by the rotation of 
both: top surface and subsurface graphene layers. Thus, 
by identifying border domains and defects, we can dis- 
criminate between moire patterns resulting from layer 
rotation on the top surface from those caused by a sub- 
surface layer rotation, and we can compare the oscillation 
amplitude of the associated line profiles. 

The experiments were performed using an STM from 
Omicron running at atmospheric pressure. The samples 
were peeled off using an adhesive tape in air. The STM 
images were collected at room temperature in both con- 
stant current and constant height mode using Ptir tips. 
The bias voltages refer to the applied sample bias, which 
corresponds to the filled-states. Every image was scanned 
from bottom-to-top. The images were flattened to cor- 
rect the sample tilting and filtered with smooth filters. 

In Fig. [ijwe show two topographic-STM images with 
atomic resolution taken on two different regions of the 
same sample: in (a) no moire pattern is observed, while 
in (b) it appears clearly. In the region without the moire 
pattern the standard HOPG image is reproduced, where 
only one lattice site is observed. The corresponding line 
profile along AB, Figure [TJc), shows an oscillation ampli- 
tude around 0.8 A. On the other hand, in the region with 
the moire pattern a triangular pattern superimposed at 
the atomic resolution is easily identifiable. This resem- 
bles giant atoms ordered in a triangular superstructure. 
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Figure 1: STM images of an HOPG sample showing the usual 
resolution (a) and a moire pattern (b). Line profiles were 
performed along AB in both cases and are shown in (c) and 
(d) respectively. 

with a corresponding line profile that has an amplitude 
of oscillation of around 1.3 A [See Figure [ijd)]. Also, in 
these oscillations we can observe the small peaks associ- 
ated to the atomic positions. 

In a given zone in the region with the moire pattern, 
we also identified a defect, a void, which provides the 
chance to observe the subsurface layer (see Fig. |2|. In- 
side this void we also observe the moire pattern, wich 
indicates that the subsurface graphene layer is proba- 
bly rotated with respect to the HOPG, while the surface 
layer is not. Additionally, other experimental evidence 
comes from comparing line profiles inside the void (LPl) 
with those and on the top surface (LP2): the line profile 
amplitude is higher inside the void. 

In order to probe the idea that the subsurface layer can 
also produce moire patterns, we calculated STM images 
for a surface model presenting a rotation in a subsurface 
layer. For this task we have performed first-principles 
total-energy electronic structure calculations [19l [20] and 
thereby obtained the spatial charge density. The elec- 
tronic density of the surface retaining its inhomogeneities 
was described by performing detailed calculations of the 
graphene system. This description is based on the ab 
initio linear muffin-tin orbital (LMTO) method, which is 
discussed elsewhere [21]. We calculate charge density on 
an energy window [Ep — e\V\^Ef] {Ej being the Fermi 
Energy) and obtained a value proportional to the current 
density detected by the STM tip when it scans the sur- 
face in constant height mode [22l [23] . To compare our 
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Figure 2: (a) Simultaneous STM observation of a void and of 
a moire pattern on the HOPG (0001) surface, (b) and (c) line 
profiles along segments LPl and LP2 highlighted on (a). 

results with STM experimental data, we introduce ra- 
dio tip effects by performing a lateral smoothing on the 
generated images. This lateral smoothing consists of a 
window average of the tunneling current density around 
a neighboring region. Our results are restricted to the 
filled states of the sample and consequently can only be 
compared with those experimental data in which the tip 
was biased positively with respect to the surface. Fur- 
ther details about these calculations can be found else- 
where [H [24] . Specifically, the model consists of a repeat- 
ing three layer graphene slab, with 14 A of separation 
between adjacent slabs to avoid electronic overlapping. 
Each slab presents a Bernal stacking sequence, and the 
central layer was rotated by 21.8°. This particular angle 
produces a commensurable superstructure [25^ composed 
of 42 carbon atoms. Figure [sj^a) shows an upper view of 
the resulting supercell while Figure |3|b) shows an upper 
view of the structure expanded by the supercell. Note 
that the top and bottom layers of our three-layer slab 
model are superimposed and the middle layer rotation 
was performed on an a-site. 

Figure [4] shows three calculated STM images for 1.0 
V bias voltage at three different values for the tip to 
sample distances: 1.0 A (left), 2.0 A (middle), and 3.0 
A (right). Brilliant (dark) regions represent high (low) 
tunneling current density. 

In all these calculated STM images we can identify a 
superstructure with a periodicity of 6.5 A. 

In addition, we observe a shifting of the density current 
maxima as was predicted before in the case of surface 



Figure 3: (Color online) (a) Ball and stick model for the up- 
per view of the supercell unit employed to perform calcula- 
tions, (b) Superstructure expanded by the supercell. Atoms 
from the top layer appear in blue, just over the atoms from 
the bottom layer. Atoms from the middle layer were painted 
green. 



layer rotation and was experimentally corroborated by 
Miller et al. [26] 




Figure 4: (Color onhne) Calculated STM images for 1.0 [V] 
bias voltage at three different tip-to-surface distances (con- 
stant height mode): 1.0, 2.0 and 3.0 A , from left to right. 
Brilliant (dark) regions correspond to high (low) tunneling 
current density. The highlighted diagonal line indicates the 
segment along which the profile in Fig. [sj^b) was obtained. 

A line profile was performed along the segment indi- 
cated on the calculated STM image for 3.0 A as tip- 
surface distance [Fig. [IJright)]. For a comparative anal- 
ysis, we have plotted such a line profile in Fig. jsjb), while 
in Fig. [sja) we have reproduced the line profile obtained 
at the same tip-to-surface distance and bias voltage for 
an array of three graphene layers with the same misori- 
entation angle, but in the surface layer [24]. Clearly, the 
oscillation amplitude associated to the surface layer ro- 
tation is larger than that associated to a subsurface layer 
rotation, thus coinciding with the experimental findings. 

Figure [6]^a) shows the band structure for the three 
layer graphene system in the usual AB stacking sequence. 
Meanwhile, Figs, ^h) and (c) correspond to the cases 
when rotations were performed on the middle and up- 
per layers, respectively. By comparing these structures 
it is possible to identify the bands origin and the effect 
of the corresponding layer rotation. Thus, the degener- 
ated band near the Fermi Level in|6ja) comes from the 
upper and bottom layers, while the single band comes 
from the middle one. So, a rotation in the middle layer 
changes the interlayer interaction and the degenerated 



Figure 5: Calculated line profiles for a three layer graphene 
array with the top layer rotated (a) and with the subsurface 
layer rotated (b). 

band looks like a monolayer graphene bandjsjb). In the 
preceding case, the middle layer contributes with a simi- 
lar band shifted in the direction of increasing energy. Of 
course a rotation in the upper layer [Fig. ^c)] lifts the 
band degeneration, and the interlayer interaction opens 
noticeable gaps among the bands. 
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Figure 6: Band structure for the three layer graphene array in 
the usual AB stacking sequence (a), and for a 21.8° rotation 
in the middle (b) and in the upper (c) layers. 

To understand the images obtained by STM, it is use- 
ful to identify each layer's contribution to the current 
detected by the STM tip. In this way, a fatband analy- 
sis for the three cases was performed and is presented in 
Figure [7| On one side, the contribution of the top atoms 
to the band structure is presented in Figure [7|^a) for the 
graphene layers in the usual AB stacking, while in (b) and 
(c) for the cases presenting rotation in the middle (sub- 
surface) and top layers, respectively. On the other side, 
the corresponding contribution of the subsurface atoms 
(middle layer) are show in (d), (e) and (f). 

The states that contribute to the STM images lie in the 
energy window [Ef^ Ef — eVbias]- With that in mind, the 
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Figure 7: Fatband analysis of the three layer graphene band 
structure. The contribution of the upper layer is represented 
for the usual AB stacking (a), for the middle layer rotation 
case (b), and for the top layer rotation case (c). The con- 
tributions of the middle layer are shown in (d), (e) and (f), 
respectively. 

following observations become interesting. Firstly, for 
the three layer graphene array without rotation [Figs. 
[7|^a) and (d)], both the upper and middle layers con- 
tribute to the STM image. There are also Van Hove Sin- 
gularities (VHS): at the Fermi Level and around 0.25 eV 
(coming from the upper layer), and near 0.5 eV (coming 
from the middle layer). This result reveals the contribu- 
tion of the subsurface layer to the STM image and gives 
new evidence to understanding the hexagonal/triangular 
structures observed by STM on the HOPG (0001) sur- 
face [6]. Secondly, a rotation performed in the middle 
layer [See Figs, ^h) and (e)] also shows contributions 
to the STM image coming from surface and subsurface 
layers. These contributions explain the superstructure 
observation, but the absence of VHS means lower tunnel- 
ing current intensity maxima and consequently, a softer 
line profile (coinciding with the experimental findings). 



Thirdly, if the surface layer has experienced the rotation, 
several VHS appear which have contributions from both 
layers. Clearly in this case there are tunneling current 
intensity maxima, which are revealed in the pronounced 
rugosity of the corresponding line profiles. Experimental 
evidence about the VHS on superstructures has appeared 
recently p7] . 

We experimentally observed the formation of moire 
patterns on HOPG surfaces. These patterns originate 
from the rotation of a graphene layer. We identified 
zones where the rotated layer is the surface layer and 
zones where the rotated layer is a subsurface one. This 
is evidenced by comparing the oscillation amplitude in 
the super-periodicity, which is lower when a subsurface 
layer rotates. The calculation of the corresponding STM 
images agrees with the experimental observation. 
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